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INTRODUCTION 


The present report gives a detailed description of a computer program for the calculation of 
the noise produced by a rotor encountering a turbulent flow field. This represents the final product 
of an extended research effort on this subject, beginning with a theoretical study of the leading edge 
noise produced by an airfoil moving through turbulence 1 . Comparison with experimental and 
further theory was given by Paterson and Amiet 2 . A more detailed examination of the effects of 
varying the various parameters in the theoretical calculation is given in reference 3. 

This computer program is critically dependent on the availability of relatively simple airfoil 
response functions that include the effects of compressibility and skewed gusts.These are given in 
references 4-8; only the necessary relations are presented here. Also needed is a turbulence model; 
that used initially was the Karman spectrum which can be found in many texts such as reference 9. 
This has more recently been extended, as discussed below. 

Th Is initial analysis was applicable to the case of the rectilinear motion of an airfoil. Using 
the principles elaborated in references 10 and 1 1, the analysis was extended in reference 12 to the 
case of rotary motion of an airfoil such as that of a propeller or helicopter rotor. This is based on an 
integration of "instantaneous spectra" around the azimuth, and includes retarded time effects on the 
generating source as discussed, for example, in reference 13. 

The analysis of reference 12 is not limited to isotropic turbulence; in principle, any 
arbitrary turbulence spectrum can be specified. Because of lack of a good non-isotropic model, 
however, the computer program was initially written assuming an isotropic turbulence. The 
present program extends this by allowing a particular type of non-isotropic turbulence. In 
particular, the turbulent spectrum input to the program is assumed to be produced by an isotropic 
turbulence that has undergone a rapid contraction. The spectrum is calculated by the analysis of 
reference 14 (together with the program noted in the following paragraph for calculation of the 
deformation tensor). This is a minor extension of the classical analysis of Ribner and Tucker 15 in 
which is calculated the spectrum of turbulence undergoing a rapid contraction. This analysis is 
based on the concept that the turbulence is composed of vorticity, behaving purely kinematically and 
moving with the local fluid velocity. The final vorticity distribution, and thus, the velocity 
spectrum, can then be determined using the classical equations of Cauchy 16 for the distortion of 
vorticity by a deformation. 

The present program does not calculate the fluid deformation that is necessary to complete the 
turbulence definition. The extension of the present program to include the non-isotropic case was 
developed in conjunction with a fluid deformation analysis of Simonich. This is a mean flow 
calculation that replaces the rotor and its wake by a series of vortex rings to simulate the flow field 
induced by the rotor tip vortices. Further details of this analysis and calculations showing the use of 
the two programs together are given in reference 17. The two programs can be used independently, 
but the deformation program, also available as an ANOPP module, is intended to be used in 
conjunction with the present noise prediction program. If this is desired, the ANOPP theory manual 
describing the deformation program (the companion to the present theory manual ), should also be 
consulted. The programs have been configured so that the output of the deformation program is 
directed as input to the present noise prediction program. 

The present program can also be run independently. For example, an isotropic turbulence can 
be assumed; the deformation tensor input is then the identity matrix (the diagonal matrix with l's 
on the diagonal). The program can also be independently run with any arbitrary deformation that 
the user wishes to specify. This allows the user to investigate the effect of the any particular 
deformation tensor on the noise. 
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ANALYSIS 


A. Airfoil In Rectilinear Motion 

The noise model assumes a vertical gust w g with the convection velocity U c parallel to the x axis. 

This gust impinges on a flat plate airfoil situated between -b < x < b and -d < y < d in the x,y plane. The 
convection velocity U c of the gust is eventually set equal to the free stream velocity U, but the initial 

derivation is performed with U c * U. Imposing the condition of no flow through the airfoil surface 

leads to an unsteady surface dipole distribution leading to radiated noise. Linearized theory is assumed 
throughout. 

The vertical gust is assumed to be produced by an incident turbulence field. The Karman 
spectrum is used as the basic model for the turbulence. This is modified later to include the effect of 
distortion produced by flow gradients. 


1. Leading Edge Noise 
A gust of the form 

« = w D (k„ky)e- i[k,(x - Uct) + k>yI (1) 

y k 

is incident on the airfoil. In general this is a single Fourier component of a more complicated velocity 
field w g (x - U c t,y) where 


w 


g (x - U c t,y) = J J # g dk*dk y 


( 2 ) 


The Inverse Fourier relation is 




JJ w g (x-U c t,»)e 


i[kx(x - U.t) + k y y] 


dx dy 


(3) 


where R is a length scale indicating the region over which the gust velocity field extends. The 
convection velocity U c can be a function of kx and k y . 


When this gust impinges on the airfoil, a surface dipole distribution is induced to oppose the gust 
flow and satisfy the condition of no flow through the airfoil surface. This surface dipole distribution 
gives a pressure jump Ap across the airfoil; a normalized pressure jump g is defined as 

Ap(x,y,t;k x ,ky) = 2 ttp 0 U c w R (k x ,k y ) g(x,k x ,k y ,M) e ^ kxUct kyy ^ (4) 


Integration over kx and k y then gives for the pressure jump due to all wavenumber components: 
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00 


Ap(x,y,t) * 2 ttPqU cj J" w Bge l(k > U ' , - |, ' <) dk.dk, 

- 00 


In the frequency domain 


00 

Ap T (x,y,co) = 2 ttp 0 J w R (X,k y ) g(x,X,k y ,M ) e -ikyy dk y 
-00 


( 5 ) 


( 6 ) 


where X is the specific value of kx such that X = q/U c and the Fourier time integral is 

T 

F T (G>)= “ 2 V / f< ° e ~ 1Wt dt T -> cd ( 7 ) 

-T 

where T = R/U c . 

Knowing the force distribution produced by the fluid on the airfoil, the far-field noise can be 
calculated by noting that a force imposed on the fluid produces a dipole pressure field response. The 
far-field pressure of a point dipole situated at x t ,yi, aligned with the z axis and with strength F z 

exp(io)t) is 


Po 


iFzO>z e i[o>t + u(Mx - ct)] e -iu[Mxi - (xxt + B z yyi)/a] 


4ttc 0 ct 2 


( 8 ) 


where the observer is located at x,y,z and 

B 2 s 1 - M 2 
u s MKx/8 2 


K x s g)/U 

a 2 5 x 2 + B 2 (y 2 + z 2 ) 


(9) 


Equation (8) gives the far-field sound produced by a point force and equation (6) gives the force 
at a point on the airfoil produced by the gust in terms of the airfoil response function. Combining and 
integrating over the airfoil planform gives the total far-field pressure. Thus, 


d b 


p( 2 >“ ) 'l^ eiI “ t * u<Mx " a)] JJ 


-iu[Mx t - (xxi + 0 2 yyi)/a] 


-d -b 


( 10 ) 


CD 


J w R (X,k y ) g(x t ,X,k y ,M) e” ,kyyi dk y dx! dy t 


-00 
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An effective lift function £. per unit span can be defined as 


b 

£(X,k y ,M) s -i-J g(x lf X,k y ,M) e -iuXl(M " x/a) dx t (11) 

-b 

This function includes retarded time effects in the propagation of the sound from the airfoil to the 
observer. For an observer overhead of the retarded source position (x/a = M ) equation (11) reduces 
to the lift/span. 

It is assumed here that there is no spanwise variation of the airfoil characteristics (or at least 
little variation on the scale of the turbulence or other upwash disturbance). The y x integration in 
equation ( 10) can then be performed. For an airfoil situated between -d < y < d (where d = semispan) 
the yi integration gives the factor 

, _ 2 sin[d(k y - uB 2 y/cr)] 

A= k _ y - liBVb? U2) 


so that from equation (10) 


P T (x,to) = i 


btozpp 
c 0 ct 2 


CO 


e i[cot + u(Mx - 



W R (X,ky) £(X,ky,M) 


A dky 


(13) 


In general, this final k y integral cannot be evaluated in closed form. Thus, the integral must be 
evaluated numerically if an exact result is desired. However, for an airfoil span large compared to the 
dimensions of the upwash disturbance a simple approximation can be made. In this case (Xd » 1 ) the 
expression in equation ( 12) varies rapidly between positive and negative values, giving a cancelling 
effect, except where 


k y = u0 z y/a = Ky (14) 

Thus, the function in equation ( 12) acts like a delta function and w 0 and £ can be moved outside the 
integral by setting k y = K y in these functions. The remaining integral of x _1 sin x is equal to tt; thus 

P T (x,co) = i ”"^ 2 Pp w R (X,K y ) «C(X,K y ,M) e ,u(Mx “ a) ( 15 ) 

where the factor exp(iut) has been omitted. This can be Fourier transformed using the inverse of 
equation (7) to find the far-field time history of the pressure; i.e.. 


CO 

p(t) = J p T (u) e lwt 

- 00 


do) 


( 16 ) 


Equations ( 15) and ( 16) give far-field expressions appropriate for the case of a deterministic 
velocity field such as the case of an airfoil cutting a vortex. The case of interest here is that of an 
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airfoil in a turbulent flow field for which the upwash velocity is not deterministic. The result for this 
case was derived in reference 1, and can be obtained from equation ( 15). Because infinite span 
airfoil response functions are used (thereby neglecting end effects) the problem is modeled as a 
turbulence region located between -d < y < d incident on an infinite span airfoil. Also, the turbulence 
is assumed to be located between -R < x < R, where R is a large length scale that can later be allowed to 
approach infinity. This avoids convergence difficulties in the integrals that would appear if R were 
immediately set equal to infinity. 

For a random function f(y) which is nonzero between -d < y < d, the spectrum function 0 is 
related to the Fourier transform of f by (see Appendix A) 

0(ky) = (TT/d) E[F d (ky) F d *(k y )j ( 17a) 

where E represents the expected value and 

d 

F d (ky) = -jL J f(y) e - ik *y dy (17b) 

-d 


Equation (15) is multiplied by its complex conjugate and the expected value taken. Then using 
equation ( 17) the following relations hold 


Thus, 


EEp T (co) p T *(<o)] = (T/tt) S pp (u) 

E[w R (X,Ky) W R *(X,Ky)j = (R d/TT*) $ ww (X,Ky) 


Spp(to) = [cobzpo/(c 0 <T z )] z ttU c d |£(X,K y ,M )| 2 $ ww (X,K y ) 


(18) 


(19) 


This agrees with equation ( 18) of reference 1. 


2. Airfoil Response Functions 

As discussed in reference 4, when calculating an approximate airfoil response function it is 
natural to divide the frequency range into a low and a high frequency regime. The ratio of airfoil chord 
to acoustic wavelength determines which regime is appropriate for a particular case. 

For simplicity the case U c = U is considered; i.e., the incident disturbance convects at the free 

stream velocity. Thus, X = Kx. The additional complexity is not needed here, but if results for the case 
U c * U are needed, they can be derived using the method of references 5 and 6. 


For low frequency and an airfoil situated between - 1 < x < 1, the normalized pressure jump is 
found to be 5 


g(x,K x ,0,M) 


_1_ f 1 - x 
ttB Y i + x 


S( K X *) e x + f(M )] 


U < 0.4 (20) 
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where 


Kx* a Kx/B 2 


f(M) s ( 1 - 8) In M + 8 ln( 1 + 8) - In 2 

Here S(k) is the Sears function 7 and the over bar on a quantity signifies normalization by the 
semichord b. 

In deriving the airfoil response function for high frequency an iteration procedure is used which 
alternately corrects the boundary condition at the leading and trailing edges. The first two terms in 
this series are 


9i(x,Kx,0,M ) 


1 e ~i©3 

tt -/tt( 1 + M )Kx( 1 + x) 


( 21 ) 


g 2 (x,Kx,0,M) » / "*■ - 1(1 + i)E*(2u (1 - x» - 1] e' i03 

tt /2tt(1 + M) Kx 


with © 3 su(l-M)(l + x) + tT/4-Kx 

and where E* is a combination of Fresnel integrals; i.e.. 


x 



0 


( 22 ) 


By comparison with numerical results of Graham® these results for g are found to be accurate to 
within a few percent if the changeover from low to high frequency is made at the value u = 0.4; the 
combination of the two functions becomes more accurate as il moves away from 0.4. 

These values are for parallel gusts; i.e., k y = 0. Values for k y * 0 can readily be found for the 
case U c - U using the similarity results of Graham®. Depending on whether the gust-airfoil 

intersection point moves subsonically or supersonically relative to the fluid, the airfoil response 
function is similar to either the 3-D incompressible solution or the 2-D compressible solution 
respectively. Results are given here for the supersonic sweep speed case, but the results for the 
subsonic case readily follow if the appropriate analytical continuations of the various functions are 
employed. 

For a gust-airfoil intersection that moves supersonically relative to the fluid, the airfoil 
response function is related to that for k y = 0 by 

g(x,Kx,K y ,M) = (8./B) g(x,Kx.,0,M„) e ]xK y Z/K * (23) 


where 

8. 2 = B 2 ( 1 + K y 2 /K x 2 ) M„ 2 = M 2 - 8 2 K y 2 /Kx 2 

Kx./8„ 2 = Kx/8 2 U. = M. Kx„/8. 2 


6 


(24) 



This relation leads to the following results for g 


g(x,K„,K y ,M) 


ttB fZ 


Ql( X,Kx,Ky,M ) * 


-x S (K«) e *[uMx + f(M.) K,*] 


-i©4 


il. <0.4 (25) 


U. > 0.4 


(26a) 


tt/ttKkU + M.)(l + x) 
gj(x,Ki,K y ,M ) - g t (x,K,,Ky,M ) / (1 ♦ x)/2 [( 1 + i) E“(2K>“ M,( 1 - x) ) - 1] (26b) 


where ©4 = Kx*[M.( 1 + x) - M 2 x - 1] + tt/4 

The effective lift Z can be calculated using equation (4) to be 

£(Kx,Ky,M) = B -1 S(Kx*) e , ** f(M * ) tJ 0 (ux/a) - iJ,(iIx/a)] 


£i(Kx,K y ,M) = 


TT Y Kx( 1 + M.) © t 


^ E*(2©i) e ’® 2 


U. < 0.4 
U. > 0.4 


(27) 

(28a) 


i© 2 


£ 2 (Kx,K y ,M) = / _ e [i(l -e" ,2@1 ) + (1 - i)[E*(4K x M.) 

IT©! / 2ttKx( 1 + M„) L 


n? 0 *i2©i 

. Y f e , , ■ E*(2u(M./M + x/ct) )] ] 

/l + (M/M„)(x/a) j 


(28b) 


where 


©1 = u(M./M - x/a) 

© 2 s X + u(M - x/a) - tt/4 

Here X = K x so that © 2 reduces to © 2 = Kx*( 1 - Mx/a) - tt/4. 


(29) 


3. Spectrum Function - Kerman spectrum 

When the forcing function is a random phenomenon such as flow turbulence certain analytic 
expressions are available to represent the spectrum. A widely used model for isotropic turbulence is 
that of von Karman. For this model the spectrum is written 


iw.(Wy,kz)= (l-k^/k J ) 


for the three-wavenumber spectrum and 


$ WW ( kx,ky ) - 


_4_ _ui Ex 2 + 

k e 2 1 + (c x 2 + lc y 2 


(30) 


(31) 
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for the two-wavenumber spectrum where 


E(k) = — 

tK } [1 + (k/k e ) 2 F /6 

I = 55 T(5/6) “u 2 

9 /W r(l/3) k | 

k _ /F T(5/6) 
e ~ L r(l/3) 


and the circumflex on fc signifies normalization by k e . Further discussion of this spectrum function 

can be found in texts such as that of Hinze*. Also, additional results pertaining to the Karman 
spectrum are given in Appendix 1 of reference 1. 
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B. Airfoil in Circular Motion 


With certain restrictions, the principles of the previous analysis for an airfoil in rectilinear 
motion can be applied to the case of a rotor. A flat plate airfoil of zero thickness can be represented by 
dipoles for both the case of rectilinear motion and circular motion (such as for a propeller). The 
far— field pressure for a dipole moving along a curved path was shown by Lowson 10 to be 


P ~ t 4tTc 0 r 2 ( 1 - M n ) 2 ~ 


1-M n 


)] 


(33) 


where the square brackets imply evaluation at the retarded time and M n s M r *n is the acceleration 

component of the dipole in the direction of the observer. The £ term represents the sound produced by 
an unsteady force and would be the same for the same force in rectilinear motion. The M n term gives 

the sound produced by a force of constant strength under acceleration. If the frequency co of the force 
is much greater than the rotational frequency Q of the rotor, the £ term is much greater than the M n 

term, and it is as if the dipole were moving in rectilinear motion at that instant of time. 

For incident turbulence noise of not too low a frequency, as well as being able to neglect the the 
M n term, the blade force amplitude for a rotor is the same as for an airfoil moving in rectilinear 

motion. The high frequencies are produced by small eddies, and small eddies are cut completely by the 
airfoil before the airfoil undergoes any appreciable rotation, and this small rotation has an 
insignificant effect on the airfoil-eddy interaction. 

Thus, the results of the previous section can be used to calculate the noise of a propeller or 
helicopter rotor at any instant of time. This result can be averaged over the period of rotation to 
obtain an average spectrum. This is a standard technique for time varying spectra as discussed by 
Bendat and Piersol 21 ; it gives a result appropriate for comparison with measurements made by a 
spectral analyzer if the measurement time is significantly greater than the blade passage time. 


The geometry of the problem is shown in Figure 1. The origin of the x,y,z coordinate system is 
fixed to the rotor hub with z along the rotor axis. The axial component of flow is in the negative z 
direction. The observer, fixed relative to the hub, is in the x,z plane at a distance r from the hub. The 
line from hub to observer makes an angle 0 with the z axis. The vector M f , representing the non-axial 

component of flow, is at an angle i|r to the y axis; for t|r = 180°, M f is along the y axis. The angle y 

represents the angle the blade span makes with the x axis; i.e., y - cot, and for y - 0 the blade span 
lies along the x axis. 

Because the blade is assumed to be a flat plate with zero steady loading, a blade segment makes an 
angle a with the x,y or azimuthal plane where 

cot a = [M t + M f cos(7 + <|f)]/M z (34) 


and M t is the azimuthal Mach number of the blade segment relative to the rotor hub. Due to the 

assumption of linearized theory, the blade segment need not be exactly aligned with the flow, but can be 
lightly loaded. The angle should not differ significantly from the value given by equation (34), 
however. 

Because of the skewed inflow, there is a spanwise component M a of flow over the blade where 



The chordwise component M b of the rotor segment Mach number is 

M b = [M t + M f cos(7 + Y)] (-1 sin 7 + J cos7) + Ic M z (36) 

The expression given by equation ( 19) for the far-field of an airfoil in rectilinear motion is 
written in terms of the present airfoil position. In order to be able to apply this relation to the 
rotating airfoil segment, the equivalent of the present position must be calculated for the airfoil 
segment. This "present position" is the position of the airfoil segment relative to the fluid if it were to 
move along a rectilinear path during the time for the sound to propagate from the source to the 
observer. An alternative method is to rewrite equation (19) in terms of the retarded source position. 
Both methods should be of equal difficulty and should give the same results. The former is used here to 
remain consistent with previous publications by the author. 

Both methods begin with a calculation of the retarded source position of the airfoil segment (or 
the rotor hub since distances on the order of the source size are not important when calculating the 
source position for a far-field observer). If at a time t = 0 a marker is placed in the fluid and a burst 
of sound is emitted, at time t = T e , when the observer hears the sound burst, the marker has moved to 

position x s . The observer is at 


x 0 = r(? sin 0 + It cos 0) 

The distance r e of the observer from the retarded source point x s is 

r e 2 = (r sin 0 - x s ) 2 + y s 2 + (r cos 0 - z $ ) 2 


(37) 


(38) 


Also, 


T e = r e /c 0 


(39) 


Since M z and Mf are given by 


M z = -M z Ic 

M f = -M f (T sin >|r + j cos i|0 (40) 


the retarded source coordinates are 


x $ = -T e c 0 Mf sin i|f = -M f r e sin Y 


y s = “Mf r e cos* 


Substitution of equation 
r e = r [ M s 


(41 ) into equation (38) leads to the following results for r e 
cos © + /I - M s 2 sin 2 © ]/( 1 - M s 2 ) 


(41) 


( 42 ) 


where 
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M s = M f + M z 


M s 2 = Wf 2 + M z 2 


( 43 ) 


M s cos 0 = Mf sin i|r sin 0 + M 2 cos 0 

By calculating (M s /M s )*(-x 0 ), © is shown to be the angle between the convection Mach number 
M s and the vector linking the observer and the source, -x 0 . Substitution of equation ( 42 ) into 
equation (41 ) determines the retarded source position x s . 

The "present” source position x p is found by adding to the retarded source vector a vector giving 
the chordwise airfoil movement during the time T e . The spanwise component of fluid velocity is 

ignored here since it was ignored in deriving equation (19). This is allowable since in calculating the 
noise of an infinite span airfoil in rectilinear motion, any spanwise flow velocity can be eliminated by 
a coordinate transformation. The only purpose of calculating the "present" source position is to find 
the appropriate coordinates for use in equation (19). This position is 

x p = x s + M.b c o T e (44) 

Introducing equations (41) and (36) into equation (44) gives 

x p /r e = -T[ M t sin 7 + M f cos7 sin (7 + ¥)] + ][M t cos7 - M f sin 7 sin (7 + i|0] (45) 


For a coordinate system at this "present" source position the observer has coordinates Xj given 
by 

Xi = Xo-x p (46) 

Equation ( 19) is appropriate for a coordinate system in which the flat plate airfoil lies in the 
x,y plane with the span along y. The Xj coordinate system above must be rotated to the same 
orientation with respect to the airfoil. Rotation of the Xj system about the z axis by an angle tt - 7 
gives a system which will be called the x 2 system with y 2 pointed along the airfoil span. Rotation of 
the x 2 system about the y 2 axis by an angle a gives the X 3 system with x 3 along the chord pointing 
from leading to trailing edge. The relations between the x t and x 3 systems are 

x 3 = (x t sin 7 - yi cos 7 ) cos a - z x sin ct 

y 3 = X!COS7 + yi sin 7 (47) 

z 3 * (x t sin 7 - yi cos 7 ) sin a + Zj cos a 

The observer coordinates in the x 3 system are the coordinates needed in equation (19). From 
equations (45) - (47) and equation (37) the observer coordinates in the x 3 system are 

x 3 = r e M t cos a - r 0 cos i 

y 3 = x 0 cos7 + M f r e sin (7 + Hr) (48) 

z 3 = (x 0 sin 7 + r e M t ) sin a + z 0 cos a 


where 


cos i = cos 0 sin a - sin 0 sin 7 cos a 
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r 0 2 = xo 2 + y 0 2 + z 0 z 


(49) 


i is the angle between the x 3 axis and the line from the observer to the rotor blade segment. This can 
be seen by replacing Xj by x 0 in the expression for x 3 given by equation (49) and substituting for x 0 
using equation (37). 


1. With No Blade-Blade Correlation 

Equations (48) - (49) give the appropriate "instantaneous observer coordinates" to use in 
calculating the sound spectrum from equation ( 19). The angle y = tot is a function of time. In order to 
calculate the time averaged spectrum, an average around the a2imuth must be determined. This 
average must account for both the Doppler shifting of frequency as the airfoil segment moves relative 
to the observer, and the retarded time effects as the blade moves around the azimuth. The latter 
correction is needed because the blade spends different amounts of time, in the acoustic sense, at each 
azimuthal location. As noted in reference 12, the proper azimuthal weighting is the factor o)/o) 0 
where 6) is the frequency of the airfoil forces and o 0 is the Doppler shifted frequency. The 
azimuthally averaged spectrum is then 

2rt 

S pp (x,co) = “ 2 “ J S' pp (x, (Oo,7 ) dy (50) 

0 


The spectrum S' pp (x,o) 0 ,7) for the case of relative motion between source and observer is related to 
S pp (x,o),7 ) with no relative motion by the Doppler factor o)/o) 0 ; i.e., as shown in reference 13, 

S' pp (x,o) 0 ,7) ■ (o)/ o) 0 ) S pp (x,o),7 ) (51) 

Thus, equation (50) and (51) together introduce a factor (o)/o) 0 ) 2 . 

The Doppler factor o)/o) 0 is shown in references 12 and 13 to be 


Uo 

0 ) 


= 1 + 


M t - SO 
1 - M r - SO 


(52) 


where M t = Mach number of source relative to observer 

M r = Mach number of source relative to fluid 

S*0 = unit vector from retarded source position to observer. 


In vector notation 


M t = M t ( -T sin 7 + J cos y ) 


M r = t ( M f sin i|i - M t sin 7 ) + j (M f cos i|i + M t cos 7 ) + lc M z (53) 

Q£= 1 (x - x s ) + ] (y - y s ) + K (z - z s ) 

Using equation (41 ) and the fact that |Q&I = r e gives the unit vector 
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(54) 


S~0 = i (x/r e + M f sin i|0 + j M f cos y + R (z/r e + M z ) 

and finally from equation (65) and S"0 = - 0 A S 

w M t [x sin 7 - M f r e cos (7 + f )] 

“So “ 1 + ( 1 - M |)r e - x M f sin i|r - 2 M z 

M t [x sin 7 - M f r„ cos (7 + ¥)] 

= 1 + — y T “ — 

/l - (M f 2 + M 2 2 ) sin 2 © 

The final result for the azimuthally averaged spectrum is given by equation (50), (51) and 
(55) with S pp (x,G),7) inequation (51) given by equation (19). The x,y,z coordinates to be used in 

equation ( 19) are x 3 ,y 3 ,z 3 , respectively, given by equation (48). 


2. With Blade-Blade Correlation 

The technique of averaging the instantaneous spectrum around the azimuth is valid if there is no 
correlation between the spectra for points widely separated in time. For the case of turbulence 
interaction, however, this is not necessarily true. The larger turbulent eddies can be cut more than 
once while passing through the rotor plane and the spectra from these events are correlated. 

For the case where blade-blade correlation exists, the autocorrelation function would consist of 
a series of peaks, such as shown in Figure 2, the n'th peak representing the correlation of sound from 
the zeroth blade with that of the n'th blade. The correlation function for the multiple blade passage 
case can be written 


Rpp(x,r,r)= 2 Rppte»*» r - nT 2 ) (56) 

n = -00 

where R pp ^ represents the cross-correlation of the far-field sound from the n'th blade with that 

from the zeroth blade. The time T 2 is the time, as heard by an observer, between blade chops of any 
given eddy. 

Figure 3 illustrates the geometry of the multiple chopping. The blades move along line AD. The 
blades are assumed to be moving rectilinearly in the vicinity of the eddy; this is a reasonable 
assumption provided that the eddies are not too large. The eddy moves in the chordwise direction 
relative to the airfoil with Mach number M b given by equation (36). After a time T t the eddy has 

moved down from the rotor plane a distance 

FE = Ti V z (57) 


and in the plane a distance 


FC = Ti V f cos(7+ 40 

At time T blade 1 is at point C; the distance BC is 


(58) 
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The distance Z is then 


BC = (T - T t ) V t 


(59) 


Z 2 = BF 2 + FE 2 = T t 2 V z 2 + [(T - T t ) V t - T t V f cos (7 + i|0] 2 (63) 

This distance is the minimum between the lines ED and BG if BE is normal to these two lines. For this 
to be true, BD 2 = BE 2 + ED 2 or 

Z 2 = T 2 V t 2 - T1 2 {V z 2 + [V t + V f cos(7 + (61) 

Equations (73) and (74) can be solved for T t and Z giving 

T t = T (V t /V z ) sin a cos a Z = TV t sina (62) 

Equation (19) for the far-field spectrum of an airfoil was obtained by multiplying equation 
( 15) by its complex conjugate and taking the expected value. The equation for the n'th blade passage is 
equation ( 15) with the additional factor exp(ik z nZ) to account for the z 3 distance nZ of the n'th blade 

path from the zeroth blade path as measured with respect to the fluid. 

The time T 2 is the time Tj plus the additional time difference for the sound from point B in 
Figure 3 to reach the observer as compared to the time for the sound from point C to reach the 
observer. From the last factor in equation ( 15), the time for the sound to reach the observer from 
point C is 


r 0 = (M b x - o)/(c o 0 b 2 ) 


(63) 


Replacing x by x + X and y by y - Y where 

XsBC = (T-T t )V t Y = T t V f sin (7 + i|0 (64) 


gives the time T t for the sound to propagate from point B. Taking the difference between T t and r 0 . 

At = X (M b - x/a)/(0 b 2 c 0 ) + Y y/(c 0 a) = T 2 - T, (65) 

Then, T 2 can be determined from equations (62), (64) and (65). 

Taking the Fourier transform of equation ( 56 ) gives 


S'pp(x» u 0 , 7) = g S'ppOi. «o,7) e ,ru0oT2 (66) 

n = -® 

The cross PSD, S pp ^ n ^ is found by multiplying equation ( 15) by its complex conjugate with the 
additional factors exp (ik z Zn) as noted before, or from equations ( 19) and (51 ) 


GO 


S'p^x.Uo.TO = J jTTU b d(a)/w 0 )|Jl(X,Ky,M)| 2 § ww (X,Ky,k z )e ik2nZ dk 2 (67) 


Introducing equation (67) into equation (66) and using the identity 

14 



£ e i2nna 
n = -os 


00 

= 25(n + a) 
n = — oo 


( 68 ) 


gives 


5 pp ( x>Gio»'y ) 


( uzpob i 


TTU b d (w/CJo) I «t(X,Ky,M )| 2 u 2 b 2 Z # ww (X,Ky,K 2 

n = -<*> 


(69) 


where 

K^ n) = (oj 0 T 2 + 2TTn)/Z 


(70) 
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C. Turbulence Contraction 


The analysis used here for calculating the spectrum of turbulence being subjected to a rapid 
contraction follows the general method of Ribner and Tucker 15 . It has been generalized slightly to 
allow for non diagonal deformation tensors, although in principle these results could be derived from 
the Ribner and Tucker results by first transforming to principle coordinates. The analysis is based on 
the fact that vorticity follows the flow. There is assumed to be no interaction of the vorticity with 
itself. The details of the analysis are presented in reference 17. A summary of the necessary 
equations is given below, and a fuller derivation is given in Appendix B. 

The upstream and downstream coordinate systems are designated by the two coordinate systems 
e u and e d respectively where e is a matrix e(i,j) composed of the three unit vectors along the axes. 
e(i,l ), e(i,2) and e( i ,3 ) denote the three components of the unit vectors parallel to the local 
vorticity vector Q, the local turbulence velocity and the unit vector in the direction of the wavevector, 
respectively. The relation between e(i,l ) u and e(j,l) d is readily obtained from the knowledge that 
vortex lines follow fluid particles. The relation is then obtained from the Cauchy relation 

=Oj d Ox f /a$p (71) 

where 3Xf/3£j represents the deformation tensor, Xj represents the upstream particle coordinates and 
£j represents the downstream particle coordinates. In the computer programs dXj/d£j is represented 

by DXDZ(I,J). The relation between e u (i,3) and e d (i,3) can be found by consideration of the manner 
of distortion of a sheet of vorticity between upstream and downstream: since vorticity follows the 
fluid, a pair of vectors lying in this plane upstream also lie in the corresponding downstream plane, 
and the wavevector is normal to the plane. First assume that e d (i,j) is given. Thus, if Qj d is replaced 

in turn by e d (i,l ) and e d ( 1,2), two vectors are found, the cross product of which lies along e u ( i ,3 ). 
The remaining vector e u (i,2) can then be found by the cross product of e u (i, 1 ) and e u (i,3). 

Knowing the relation between e u (i,j) and e d (1,j), the relation between the downstream and 
upstream wavenumbers is given in Reference 17 as 

k d /k u = e u (i,3) e d (j,3) 3Xj/35j (72) 

The relation between fluid velocities is 

qj U (k u ) = q d (k d )f 1 u e(i,2) k d /k u (73) 

where fi u is the magnitude of the vector found by replacing Qj d in equation (7 1 ) by e( j,l ). 
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List of Symbols 


B 

b 

c 

c 0 

d 

E[.J 

E(k) 

E* 

u 

f(M) 

9 
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1. 1, K 
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k x , k y 
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Kx 

K y 
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M 

M a 

-f 

Mr 

-t 
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Pt 

Ap 

Po 

Q 

fl 


r 0 


r 



(n) 



T 

Ti 


Blade number 

Airfoil semichord 

Airfoil chord 

Sound speed 

Airfoil semi span 

Expected value operator 

Energy spectrum; see Eq. (32) 

Fresnel integrals; see Eq. (22) 

Frequency measured by observer 
Function of Mach number; see Eq. (20) 

Normalized airfoil response function; see Eq. (4) 

Factor defined by Eq. (32) 

Unit vectors in the x, y and z directions 
Bessel functions 
x and y wavenumbers 

Wavenumber defined by Eq. (32) 

U/U 

K x /B 2 

uB 2 y/a 

Turbulence integral length scale 
Effective lift; see Eq. ( 1 1 ) 

Free stream Mach number 
Spanwise flow Mach number over blade 

Chordwise Mach number of blade 
Flight Mach number in rotor plane 
Mach number of rotor segment relative to fluid 
Mach number of blade segment relative to rotor hub 
Flow Mach number along rotor axis 

Unit vector from retarded source to observer 
Fourier time transform of p 

Pressure jump on airfoil produced by all gusts 

Far-field pressure of dipole 

Magnitude of turbulent velocity in Appendix B 

Turbulent velocity in Appendix B 

Distance from observer to retarded source 

Distance of observer from rotor hub 
-x 

Integration limit for spatial Fourier transform 
Autocorrelation of far-field sound 

Cross correlation of far-field sound from n'th blade with that of zeroth blade 
Cross spectrum of airfoil surface pressure 

Far-field acoustic spectrum for observer fixed relative to blade segment 
Far-field acoustic spectrum for observer fixed relative to blade hub 

Time between blade passes 
Time between eddy chops 
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T e 

t 

u 2 
U 
Ub 
u„ 


w g 

W R 

x,y,z 

x 0 
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*p 

Xi 

X3 

X 

y 

z 

a 

8 2 

Bb 2 

r 

7 

© 
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6 

X 


U 

Po 

a 2 

r 0 

$ 

*ww 

V 

(0 


0)q 

(0 

*•»» 

Overbar 

Subscript 


Time between eddy chops as heard by observer 
Time between emission and reception of sound 

Time 

Mean square turbulence level 
Stream velocity 

Component of fluid velocity along airfoil chord 
Turbulence convection velocity 
Flow velocity in rotor plane 
Azimuthal velocity of airfoil segment 
Axial flow velocity 
Gust velocity 

Spatial Fourier transform of gust velocity 

Cartesian coordinate system; flow along x and y along span 

Vector from rotor hub to observer 

Vector from rotor hub to retarded source position 

Hypothetical "present source position" defined by Eq. (44) 

Vector from "present source position" to observer 
Observer position in rotated coordinate system 
Distance defined by Eq. (64) 

Distance defined by Eq. (64) 

Distance defined by Eq. (62) 

Out of plane angle of rotor blade 

1 - M 2 

1-M b 2 

Gamma function 
Azimuthal blade angle 
Angle defined by Eq. (43) 

Functions defined by Eq. (29) 

Angle of observer from upstream rotor axis 
to /U c ; also, in Appendix B this is the wavelength 

MK x /8 2 

Density of air 
x 2 + 0 2 (y 2 + z 2 ) 

Time between eddy chop and when observer hears it 
Angle defined by Eq. (49) 

Velocity spectrum of turbulence 

Angle of flow in rotor plane; see Fig. 1 

Radian frequency measured by observer fixed relative to rotor hub; also, vorticity in 
Appendix B. 

Doppler shifted radian frequency measured by observer fixed to rotor hub 
0)5 /U 

Fourier transformed variable 

Normalized by semichord b 

Denotes similarity variable; see Eq. (28) 
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RppendiH ft: Relation of Fourier Transform of a Random Function to Its 

Spectrum 

If f(y) is a random function of y between -« < y < «, the infinite Fourier transform between 
these limits is not defined, in general, because of convergence problems. If, however, finite limits are 
placed such that f ( y ) is random between -d < y < d and f( y ) is zero outside this range, a Fourier 
transform 

d 

F d (k y ) = -^- J e-M f(y) dy (Al) 

-d 

can be defined for any particular sample f(y). The cross correlation function R is 

R(il) = E[f(y + n) f(y)l (A2) 

The spectrum 0 is the Fourier transform of R so that 

00 00 

0(k y ) = R(ti) e _ik V 11 d H= 27 j E[f(y + n) f(y)l e" ik ^ dn 

- 00 - oo 

d 

= ~2r\ E [ f(y) J f( C )e ’ k(y " C) dC ] (A3) 

-d 

= E [F d (k y )e ik yy f(y) ] 

Integration over y then gives 

2d 0(k y ) = 2tt E [ F d (k y ) F d *(k y ) ] (A4) 

which is equation ( 17). 
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RppendiH B: The Spectrum of Turbulence Undergoing a Rapid Distortion 


The following analysis calculates the change in wavevector and fluid velocity for a turbulent flow 
undergoing a rapid contraction. A sketch of the problem is shown in figure 4. The analysis begins in 
Section 1 by deriving the relations between vorticity and velocity for a particular wavevector 
component at a point. These results are derived from the assumed sinusoidal variation for the velocity 
field of a wavevector component. 

In section 2 the equations for the transport of vorticity are used to derive a relation between the 
pre-contraction and post-contraction velocity vectors of a wavevector component. In section 3 the 
same transport equations are used to derive a relation between the pre-contraction and 
post-contraction wavevectors. 

Finally, in section 4 the results are compared to the results of Ribner and Tucker. Agreement is 
found for both the wavevector and the fluid velocity transformation. 


1. Relations Between Velocity and Vorticity 

In calculating the noise generated by a rotor moving through a turbulent field, a required input 
is the turbulence spectrum as a function of the wavevector 

k = i k x + j k y + K k z (Bl) 

The velocity field g. of a single wavevector Fourier component can be written 

g d (k d , x) = Q d (k d )e- d ‘- (B2) 

where the d superscript refers to the downstream or post-contraction location, at the rotor face. The 
superscript u will refer to the upstream pre-contraction velocity field. The object of this section is 
to relate the two velocity fields. In equation (B2) Q and k are orthogonal from the assumption of 
incompressible flow (see figure 5). Thus, Q can be written as the cross product of k with some vector 
P; i.e.. 


Q(k) = k x p(k) 


( B3 ) 


The superscripts u and d are not used here since the equation can be applied to either region. 
The vorticity field g) of the wavevector Fourier component is given by 


6) = V*3 


(B4) 


so that from equations (B2) and (B3) 


gj(k,x) = ik x Q(k ) e - * - = i[k(k • P ) - P k 2 ] e - * - 


(B5) 


From equation (B5) it will be noted that 


21 



Writing 


k • to = 0 


(B6) 


co(k,x) = Q(k) e’- * - (B7) 

Q is to be solved for in terms of Q. Then, knowing how Q transforms between pre-contraction and 
post-contraction locations, the behavior of the velocity field will be known. From equations (B5) and 
(B7) 

Q(k) = i [k(k • P ) - Pk 2 ] (B8) 

On taking the cross product of this equation with k, the first term on the right hand side drops out 
leaving 

k * Q = -ik 2 k * P = -ik 2 Q (B9) 

This gives two equivalent expressions for relating vorticity and velocity of a Fourier component 

G=ik*G Q = i kxQ/k 2 (B10) 


2. Relating Pre-contraction and Post-contraction Velocities 

The relations between pre-contraction and post-contraction are now considered. The vorticity 
of a fluid particle is related by the deformation tensor. Using an expression attributed to Cauchy 16 

* Qj u d5j/dXj Qj u = Qj d ax/aSj (Bll) 

where x { represents the pre-contraction coordinate system and the post-contraction coordinate 

system. Both coordinate systems are assumed right-handed and orthogonal. The coordinates of a fluid 
particle also transform according to equation (BID. For example, consider a cube of 
post-contraction fluid with edges parallel to the axes. The edges can then be written (11,0,0 ) , (0,8,0), 
(0,0,8) where 8 is the length of the side of the cube. By writing each of these vectors as a row vector, 
the three together form a 3*3 matrix C d which can be written 

C d = 8 1 (B 12) 

where I is the unity matrix. The three upstream transformed vectors denoting the cube are then given 
by 

c u = c d [ax i /ae j ] =8[ax i /as j ] (Bi3) 

where [3x^/3^] denotes the deformation tensor defined as 
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axf/asj = 


(B 14) 


axj/a^i 3x1/352 3x1/353 

3x2/35 1 3x2/352 3 X 2/35 3 

3x3/351 3x3/352 3x3/353 

The pre-contraction volume of the cube is then given by the cross product of two of the vectors with 
the dot product of the third, which is the same as the value of the determinant of the matrix. For 
incompressible flow this must be J 3 , the same as the post-contraction value. This requires that the 
determinant of the deformation matrix in equation (B14) be 1. 

The following paragraphs will give the detailed procedure for transforming between 
pre-contraction and post-contraction velocity fields. One basic assumption of the analysis to follow is 
that fluid planes in one coordinate system remain planes after deformation. This assumption will be 
grossly incorrect, in general, on a large scale. However, it becomes more and more accurate as 
restriction to smaller and smaller scales is made. The assumption should thus be adequate if the 
turbulence scale is small compared to the scale of the distorted flow. The assumption arises when a 
one-to-one correspondence is made between a pre-contraction and a post-contraction Fourier 
component of the turbulence. Each Fourier component consists of a vorticity distribution over all 
space, whereas the deformation tensor is a function of position, and the flow is distorted by different 
amounts at different points. 

The relation between the pre-contraction and post-contraction wavevectors has yet to be 
determined. For the present they will be denoted by k u and k d with the assumption that the two 
wavevectors are related by the deformation tensor in some manner to be determined later. Equations 
(BIO) and (Bll) can then be combined to obtain a relation between the pre-contraction and 
post-contraction velocity field. Equations (BID relate pre-contraction and post-contraction 
vorticity levels while equations (BIO) relate velocity to vorticity at either pre-contraction and 
post-contraction location. 

Combining equations (BlOb) and (BDb) gives 

Q i u (k u ) = i (kj u /k 2U ) Q n d (k d ) (3x k /35 n ) £ jjk (B 15) 

where Cartesian tensor notation is used. Summation over repeated indices is assumed, and €jj k is the 
alternating tensor. Using equation (B 10a) to substitute for Q n d gives 

k,'V(k'i) £ 1mn Ox k /35 n ) £ ijk (B 16) 

This equation relates the pre-contraction and post-contraction velocities if the relation between the 
pre-contraction and post-contraction wavevectors k u and k d is known; this relation will be 
determined shortly. It should be noted that although the spatial coordinates of the pre-contraction and 
post-contraction locations do not appear explicitly, they do appear implicitly through the deformation 
tensor. Thus, even though the velocity field has been decomposed into spatial Fourier components, this 
is in some sense a local decomposition. 

In applying equation (B 16), in addition to specifying the deformation tensor and wavevector of 
interest, the velocity Q d must be specified. The amplitude of Q d is not of concern since the ratio of 
amplitudes is all that is needed. The direction of Q d must be specified, however. Once k d is specified, 
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Q d can lie anywhere in the plane normal to k d since k d and Q d are normal for the velocity modes 
considered. For the specific case of an airfoil moving through the turbulent field consider a plane 
formed by the vector k d and a unit vector n normal to the airfoil surface. Then for the specific k d 
vector chosen, the velocity field can be decomposed into a component in the plane formed by the vectors 
k d and n and a component normal to this plane. In calculating the airfoil response to turbulence using a 
linearized analysis, only the velocity component along n will give a contribution. Thus, any component 
of O d normal to the k d , n plane can be neglected. This allows Q d to be specified as normal to k d and in 
the k d , n plane. 

3. Relating Pre-contraction and Post-contraction Wavevectors 

The only remaining relation to be determined in equation (B 16) is that between k u and k d ; i.e., 
how the wavevectors of the mode considered become distorted in going from upstream to downstream. 

As with Q u and Q d , the wavevectors will be related through the deformation tensor. The three vectors 
Q d , Q d and k d form an orthogonal system. In addition, k d is assumed to be in the direction Q d * Q d so 
that the system is right handed. Define vectors ii d » u , e 2 d » u , e 3 d,u along Q d * u , O d » u , and k d » u 
respectively where either d or u is read, not both. First the relations between e u and e" will be 
determined and then the relation between k u and k d . As a matter of notation, note that subscripts here 
can refer to either components of a vector as in the Cartesian tensor notation of equation (B 1 6 ), or to 
a specific vector such as e t , e 2 , e 3 . Generally the meaning will be clear from the context since e 
already indicates a vector and the subscripts 1,2 or 3 then indicate a specific vector, not the 
component of a vector. 

The relation between ej u , e t d is readily obtained from the knowledge that vortex lines follow 
fluid particles. Since e t is parallel to Q, equation (Bll) immediately gives e t d ; i.e., if fj defines the 
direction of e^, then 

(fi u )j = (ei d )j (dx/a^j) e t u =_f, u /|ii u | (B 17) 

where the j subscript on e d now indicates one of the three Cartesian components of the vector £t d . 

Also, the summation convention over repeated indices is assumed. 

To calculate the vector e 2 u , it is not sufficient to merely substitute subscripts 2 for 1 in 
equation (B 17). The resulting vector so defined would not necessarily be orthogonal to ej u . Therefore 
a different approach will be used, first calculating e 3 u , then e 2 u . First define a vector 

(f 2 u )j = (e 2 d )j Ox/aep (B 18) 

As noted above, this vector will not in general be orthogonal to ej u . However, it will lie in the same 
plane of vorticity as ej d and e 2 d . For example, consider a plane of vortex lines defined by the vectors 
e t d and e 2 d . Because the vorticity moves with the fluid and because equation (B 18) is a Lagrangian 
type of equation following fluid particles, f 2 u must lie in the same plane of vortex lines. Thus, 
although e 2 u hasn't yet been determined, e 3 u can be found since it must be normal to the plane defined 
by e t u and f 2 u ; i.e., 


eV = eV*f 2 u /|e 1 Ux _f 2 u | (B19) 

where, as in equation (B 17), the denominator is for the purpose of normalizing the result to a unit 
vector. Finally, e 2 u can be found directly from e t u and e 3 u ; i.e.. 
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e* u -e 8 u «ei u (B20) 

Now the relation between k u and k d can be determined. This relation is found from the relations 
between wavenumber k and wavelength X 

k u = 2ir/X u k d =2tT/X d (B21) 

along with the value ofi 3 U obtained by replacing k d with Qj d in equation (B 1 la) and by the knowledge 

that k u and k d are parallel to the vectors e 3 U and e 3 d respectively. Consider a vector X d e 3 d in the 
post-contraction fluid. If this vector begins on the crest of a wave of the turbulence wavevector 
component, it will reach just to the next crest since the vector has magnitude X d; see figure 6. In the 
pre-contraction fluid the vector begins and ends on the corresponding adjacent crests, but not 
necessarily along the shortest distance; i.e., the vector will not necessarily be along e 3 u . If this 
upstream vector is denoted h 3 U , then following the definitions in equations (B17) and (B18) 

(h 3 u )j = X d (f 3 u )j = X d (e 3 d )j Ox/^j) (B22) 

The component of h 3 u in the direction of e 3 u must have a length X u from the above discussion. Thus, 

X u = h 3 u . e 3 u = X d (e 3 u ) i (e 3 d )j (dXj/d$j) (B23) 


From equation (B21) 


k d /k u = (e 3 u )j (e 3 d )j Ox/aSj) (B24) 

The relations between the pre-contraction and post-contraction values of both wavenumber and 
velocity can be determined from the above relations. The procedure is first to find the three 
pre-contraction e u vectors from the three post-contraction e d vectors using equations (B 17) through 
(B20). The ratio k d /k u can then be found from equation (B24). Finally, Q u can be found from 
equation (B 16). 

Several of the vector operations in equation (B 16) need not actually be carried out. The 
necessary operations have already been performed in deriving the e u vectors. Thus, since kj d and Q m d 

are orthogonal. 


k l' 1 °mNmn - - < e . d >n kd 0,1 < B 25> 

The vector (e^j dx t /d£j has already been calculated in equation (B 17) and is the vector _fi U parallel 
to e t u . Thus equation (B 16) can be written 

O i u (k u ) = Q d (k d ) f t u e 2 u k d / k u (B26) 

A computer program was initially written using equation (B16). This was changed to the simplified 
version using equation (B26), and the programs were found to agree. 


4. Comparison with Results of Ribner and Tucker 
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a. Wavenumber Relation 


Ribner and Tucker 15 have performed this same analysis, but for the restricted case where the 
deformation tensor is diagonal. This is no restriction in principle since a rotation to principal 
coordinates can always be performed in which the deformation tensor is diagonal. In practice the 
present method may be simpler. The present results will be examined to verify that they reduce to 
those of Ribner and Tucker. 

The deformation considered by Ribner and Tucker is 


r ax "1 

nr 1 

0 

0 


II 

E 

0 

_0 

nr 1 

0 

0 

is* 1 

(B27) 


The relation given between k d and k u is 

k d = ( ki u /l,, k 2 u /fi 2 , ks u /«3 ) ( B 28 ) 

This can be shown to follow from equations (B24) and (B27); i.e., multiplication of equation (B24) 
by k d gives 

k 2d = k^ kj d ax i /ae j = k t u k t d Jti" 1 + k 2 u k 2 d flz" 1 + k s u k 3 d (B29) 

This will hold for all k values if 

ki u ar 1= k t d k2 U fl2- 1 = k 2 d k 3 u J 3 " 1 = k 3 d (B30) 


Equation (B30) agree with equation (B28). Thus, the procedure for calculation of k u from k d 
presented here reduces to that of Ribner and Tucker for a diagonal deformation tensor. 


b. Velocity Relation 

Equation (B 16) for the relation between pre-contraction and post-contraction velocities can 
similarly be reduced to the form given by Ribner and Tucker. Notation presents a minor problem in 
that the use of the Cartesian indices can be confusing. The notation used here is that a subscript in 
parentheses is a duplicate index and is not summed over unless the index appears twice elsewhere in 
the expression. Thus, A i B/q implies no summation over i, but the variable B has the same index as A. 

However, for Bj C( 1 ), summation over i is assumed with C in each term of the summation taking the 
same index as A and B. 

Equation (B 16) is rewritten into the form used by Ribner and Tucker relating the 
post-contraction velocity Q d to the pre-contraction value Q u (rather than vice versa as in equation 
(B 16) ). Then 


0,«<k<i) - - (yf/k*) k,"0 m “(k") £ , mn 0| k /x„) £|)k 


(B31) 


Using the inverse of equation (B27) gives 
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Q i d (k d ) = - (kj d /k 2d ) k l u Q m u (k u ) £ 1mn c ijn l (n) 


(B32) 


Now 

ii 82 = 1 

(B33) 

for incompressible flow since the determinant of the deformation matrix must be 1 as discussed 
previously. Then since Cjj n = 0 unless i * j * n, equation (B32) can be written 

O, '•<£"> - - U/Vk* 1 ) k, u Q m u (£“) « 1mn e i j n H(i) «<)> l ’ 1 
Using the relation 

(B34) 

e lmn c ijn = ®il ^jm ” ^im ®jl 

(B35) 

equation (B34) becomes 


Q 1 d = (k j d /k 2 d )[ 8 (i) l (j) ]'* (kj u Q , u - k i u Q j u ) 

(B36) 

Introducing equation (B30) 


u .Ufc u 

q,i. ,* ro, u - 1 1 d o, u i 

1 ‘(o L Mi) k J 

(B37) 

But this is exactly equation ( 13) of reference 15. Thus, the present result relating 
and post-contraction velocity reduces to the previously derived result of Ribner and 
case of a diagonal deformation tensor. 

pre-contraction 
Tucker for the 
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Figure 2: Samp!* cross-correlation of far-fWd note*. 
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Hgure 3: Kinematic* of eddy cuUtoQ pwcaaa. 


30 


i 

I 





32 
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e 3 d transforms to a vedor fi using equation (B1 1). 
X u is found from the relation X u -B* e 3 u 


Figure 6: Relation between pre-contraction and post-contraction wavevectors. 



Report Documentation Page 

Space A<ynrtSra*on 

1. Report No. 

NASA CR-181788 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Tide and Subtitle 

Noise Produced by Turbulent Flow Into a Rotor: 
Theory Manual for Noise Calculation 

5. Report Oate 
June 1989 

6. Performing Organization Code 

7. Author(t) 

R. K. Amiet 

8. Performing Organization Report No. 

10. Work Unit No. 

505-63-51 

! 

9. Performing Organization Name and Address 
United Technologies Research Center 
East Hartford, CT, 06108 

11. Contract or Grant No. 
NASl-17763 


13. Type of Report and Period Covered 
Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 

14. Sponsoring Apancy Code 

15. Supplementary Notes 1 


Langley Technical Monitor: 
Final Report 


Thomas F. Brooks 


16. Abstract 

An analysis is presented lor the calculation of noise produced by turbulent flow into 
a helicopter rotor. The method is based on the analysis of Amiet for the sound pro- 
duced by an airfoil moving in rectilinear motion through a turbulent flow field. The 
rectilinear motion results are used in a quasi-steady manner to calculate the instan- 
taneous spectrum of the rotor noise at any given rotor position; the overall spec- 
trum is then found by averaging the instantaneous spectrum over all rotor azimuth an- 
gles. Account is taken of the fact that the rotor spends different amounts of re- 
tarded time at different rotor positions. Blade to blade correlation is included in 
the analysis, leading to harmonics of blade passing frequency. The spectrum of the 
turbulence entering the rotor is calculated by applying rapid distortion theory to 
an isotropic turbulence spectrum, assuming that the turbulence is stretched as it is 
pulled into the rotor. The inputs to the program are obtained from the atmospheric 
turbulence model and mean flow distortion calculation, described in another (Volume 
of this set of reports. This study provides the analytical basis for a module^which 
has been incorporated in NASA’s ROTONET helicopter noise prediction program. 


17. Key Words (Suggested by Author(s). 

turbulence noise, 
rotor noise, nar row band noise 


18. Distribution Statement 
Unclassified - Unlimited 


Subject Category 71 


19. Security Ctassif. (of this reportl 

20. Security Ctassif. (of this page) 

21. No. of pages 

22. Price 

Unclassified 

Unclassified 

35 

A0 3 


NASA FORM 1S20 OCT 86 






















